Topological transport in a spin-orbit coupled bosonic Mott insulator 
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We investigate topological transport in a spin-orbit coupled bosonic Mott insulator. We show 
that interactions can lead to anomalous quasi-particle dynamics even when the spin-orbit coupling is 
abelian. To illustrate the latter, we consider the spin-orbit coupling realized in the experiment of Lin 
et al. [Nature (London) 471, 83 (2011)]. For this spin-orbit coupling, we compute the quasiparticle 
dispersions and spectral weights, the interaction-induced momentum space Berry curvature, and the 
momentum space distribution of spin density, and propose experimental signatures. Furthermore, 
we find that in our approximation for the single-particle propagator, the ground state can in principle 
support an integer Hall conductivity if the sum of the Chern numbers of the hole bands is nonzero. 



Introduction-Spm-orbit coupled electronic systems 
that break time reversal symmetry can exhibit the quan- 
tum anomalous Hall (QAH) phase, which supports an 
integer Hall conductivity in the absence of a magnetic 
field, and was first proposed in a model for grapheneP 
At the single particle level, this is due to the nontrivial 
momentum space topology of the Hamiltonian, which is 
characterized by a momentum space Berry curvature that 
causes in the semiclassical equations of motion an anoma- 
lous velocity transverse to external forces, and its integral 
over the Brillouin zone is a topological invariant known as 
the Chern number. 2 Since it is now possible to engineer 
spin-orbit couplings in cold atomic systems P there have 
been several recent proposals of cold atomic optical lat- 
tice realizations of the QAHphasep^and the related atom 
topological insulators (TI), 5 which can be built from two 
time-reversal conjugate QAH systems P^All these pro- 
posals concern fermions. 

While TTs and QAH phases have been extensively 
studied at the non-interacting single particle level and 
in the presence of weak interactions, the strongly inter- 
acting regime is much less understood. This regime is 
relevant for some strongly correlated materials such as 
the transition metal oxides,^ and can be accessed in cold 
atomic optical lattices P which has the advantage that one 
can study interaction effects in a perfectly clean environ- 
ment and create topological phases for bosons, which we 
will study in this paper. 

So far, only abelian spin-orbit couplings have been 
achieved in cold atoms, as in the experiments of Ref. [3]. 
This spin-orbit coupling, equivalent to an equal amount 
of Rashba and Dresselhaus coupling, has zero Berry cur- 
vature in the absence of interactions. In this paper, we 
consider the effect of an optical lattice and interactions 
in the Mott-insulating regime for this system. We find 
that the interactions renormalize the momentum space 
spin texture and generate Berry curvature, resulting in 
an anomalous velocity for the quasiparticles. 

In two dimensional topological insulators such as HgTe 
quantum wells, one can change the topological phase by 
tuning a time-reversal breaking gap in the single parti- 
cle Hamiltonian!^ In a spinful Mott insulator one may 
consider whether such a gap can arise due to interac- 



tions. We find that this is indeed possible, and that the 
Chern number can change in the (ferromagnetic) Mott- 
insulating phase. If the sum of the hole band Chern num- 
bers is nonzero, we find, in the approximation that the 
propagator contains only quasiparticle poles, the ground 
state can support an integer Hall conductivity. 

Interacting Berry curvature and Hall conductivity- For 
an interacting system, the Hall conductivity an can be 
expressed in terms of the momentum space single par- 
ticle propagator G a ^(o;,k), wh ich r eads in two spatial 
dimensions at zero temperature^^ 
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where here and below i = (k x , k y ) and + denotes a pos- 
itive infinitesimal, h is the Planck constant (this con- 
ductivity gives a mass current as the atoms are neu- 
tral). This expression is a topological invariant in fre- 
quency and momentum space^ called the Chern number 
for an interacting system and is related to the number 
of edge states, and a nonzero integer value defines QAH 
phaseP^E3l Next we show that this phase can exist for a 
bosonic Mott insulator that supports quasi-hole excita- 
tions with non-trivial Berry curvature, as is typically the 
case in our approximation. 

Deep in the Mott insulating phase, the relevant excita- 
tions are long-lived quasiparticle/quasihole states whose 
spin-orbit coupling is encoded in the Hermitian matrix 
structure of G, which can be expressed the local spin ba- 
sis as G(cj,k) = UGdU^, where f7(cj,k) is a unitary ma- 
trix and Gd is a diagonal matrix. Below we will compute 
G(o;,k) in the Mott insulator. Using this basis, picking 
up the contributions from the quasiparticle poles in Gd in 
the uj integration in Eq. 0, and making the finite tem- 
perature generalization of Eq. 0, the Hall conductivity 
can be expressed as o~h = ^2 sn v sn /h, where 
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where T is the temperature, ks the boltzmann con- 
stant, s, n = zb are indices for the spin-orbit and the 
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quasi-particle/hole bands, respectively, and the interact- 
ing Berry curvature is given by 15 
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where v ksn = <9 k ek sn are the band velocities, and 6ksn &re 
the quasiparticle energies relative to the chemical poten- 
tial. The Berry electric field £ and magnetic field 
which in 2D points in-plane and out-of-plane, respec- 
tively, can be expressed in terms of a field strength tensor 
{Bl,£l) = (J^V-F^), where T% = - d v A^ with 

/i, v = (uj,k x ,k y ) is defined in terms the diagonal com- 
ponents of the matrix gauge field A^ = [iU^d^U] ss . It is 
convenient to define on-shell matrix rotations and gauge 
fields in k space by U\^ sn = U(uj = ek sn ? k) and 

A^sn = i[U ] dJj} ss = [,4>ksn + >A k ]|u;= eksn , 

then with analogous definitions as above, B^. sn = T^ y . 

In the zero temperature limit, only the quasihole con- 
tributions in Eq. Q with ek sn < remain, 
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where the integral over the Brillouin zone (BZ) is the 
Chern number, which is a topological property of each 
band. If the sum of these numbers is nonzero, this results 
in a ground state, integer Hall conductivity. Thus, al- 
though the QAH phase cannot occur for non-interacting 
bosons, it can occur in the bosonic Mott insulator, or 
more generally, for any bosonic system with hole bands 
of which the sum of Chern numbers is nonzero. 

Spin-orbit coupled Bose Hubbard model- In the remain- 
der of this paper we compute the Green's function in the 
Mott insulator and we specifically consider a two dimen- 
sional (2D) square optical lattice of bosons with two spin 
components deep in the Mott-insulating phase at a com- 
mensurate filling, with spin-dependent hopping ampli- 
tudes. We will describe the system by a pseudo-spin 1/2 
Bose Hubbard model with onsite repulsive interactions, 
and write the Hamiltonian as H = Hq + V, where 



H = U a(3<4a a lf3 a if3 a ia 5 
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where i denotes lattices sites, <> denotes nearest neigh- 
bors, ce, /3 = zb are spin indices, hats denote matrices in 
spin space and a k and ai denote two-component spinor 
field operators in momentum and real space, respectively, 
and tij are hopping matrices. It is convenient to write 
^k = ^k + h k • 0", where & is the vector of Pauli matri- 
ces, and h k is the non-interacting spin-orbit field. The 
dispersions of the non-interacting particles are given by 
eika = ftk + a|hk|. We will first derive results for a gen- 
eral h k , and then apply them to the spin-orbit coupling 
of Ref. [3]. 
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FIG. 1: For the filling fraction (1,0): Dispersions with energy 
in units of U with n = A = .01 for (a) spin up and down quasi- 
particles and (b) spin up quasiholes. (c) Ground state momen- 
tum distribution of spin density Sk for rjy/m/U = A/U = .04. 
The vector field represent the in-plane components and the 
colored density plot represent the z-component. (d) The 
ground state Berry electric field 6 X and interacting Berry cur- 
vature B z of the spin up quasiparticle band, plotted as func- 
tion of dimensionless wave vectors in the region \k' X)y \ < 0.2. 



Strong coupling perturbation theory-We will compute 
the single particle green function in the strong coupling 
limit V/U a p <C 1, using a perturbation theory in which 
the hopping Hamiltonian V is the perturbation to the in- 
teraction Hamiltonian ff Q ! 16 l 17 l To this end, consider the 
grand canonical partition function with external sources 
J, which has the path integral representation (here and 
below we set h = 1) 
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where S = f dr Y,ial a ia( d r ~ + H (a)] and r 

denotes imaginary time. We will use it as a gener- 
ating functional for the correlation functions given by 

(ai a a]p) = (l/Z)(5 2 Z/5J %a 5J*p)\j =0 . 

We will specify the unperturbed eigenstates in the oc- 
cupation number basis by the number of spin up and spin 
down particles per site, (iY+,iV_) with energies per site 
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given by 



U a 



N a (N a - l)-fi a N a 
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The unperturbed imaginary-time ordered corre- 
lation function is defined by SijSa(3g(3(T,T f ) = 
-(N+N-\Ta ia (T')(J jl3 (T)\N+N-), which is local in 
space and diagonal in spin, and its Fourier transform 
reads 
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where £q,(7V + , 7V_) = £"iv a +i,iv_ a — EN a ,N- a are the ex- 
citation energies, and uj n are bosonic matsubara frequen- 
cies. We will only be interested in the zero temperature 
limit, and the dependence on N a above has been cal- 
culated in this limit. We also define the total particle 
go and spin g z components by g = (g + +g-)/2,g z = 

(g+-g-)/2. 

Carrying out this perturbation theory, we find 

- G-^iwn.k; JV a ) = -g-\iu n ;N a ) + h(k) . (9) 

Although this expression is first order in hopping, its in- 
verse, the propagator, contains all orders in hopping. 

In the following, we will denote by the replacement 
iuj n — » uj the analytic continuation of the green func- 
tions, which will be necessary to compute the spectral 
functions and quasiparticle dispersions, and to define the 
spin texture in (cj,k) space. Defining components by 
— G _1 = d + d • 6", the spin textures are given by 

d(oj, k; N a ) = -[g-\uj- N a )] z Z + h(k) , (10) 

and d(w, k; N a ) = — [g~ l (uj] N a )]o + h(k). The in-plane 
texture is unchanged, while the z component is shifted 
by the a z component of the inverse propagator. 

The quasiparticle dispersions are computed by solving 
— (G~ 1 1 ) ss (uj = ek S n) = d + s|d| = 0. The Berry electro- 
magnetic fields are given by 



s - 



. (11) 



where d = d/|d| and the expressions above are evaluated 
at uj = eksn- Defining the quasiparticle spin-orbit field in 
k space by fksn = d(uj = ek S n,k), it is readily verified 
that Sk sn /27r = 0/47r)f ksn • dj^ sn x ^f ksn , and thus 
the Chern number is the integer winding number of fksn- 
From the simple poles of G = (d — d • cr)/(d 2 - d 2 ), 
we find that quasiparticle peaks in the spectral function 
A(uj) = i(G(uj + i0 + ) - G(uj - z0 + )) are given by, 



i(cj,k) = -27T^Z ks 
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where Z\^ sn = [d u (d + s\d\)\ u = e ^ sn ] 1 are the quasipar- 
ticle weights. The numerator above is a spin projection 
operator that can be written as XksnxLm> wnere Xksn 
is the normalized spinor part of the quasiparticle wave 
functions which points in the direction fksn- We have 
checked that the trace of the spectral function contains 
only these quasiparticle peaks. The momentum distribu- 
tions of particle and spin density are given by {rik, Sk} = 
f(duj/27r)tr[{l,&}A(uj,k)}/(e UJ / kBT - 1), which in the 
ground state (cf. Eq. Q) contains only quasihole contri- 
butions given by {n k ,s k } = E s ^ks-{l, sfks-}- 

Next, we consider when the interaction induced spin- 
orbit texture fksn can have a nonzero winding number. 
For this to be true, hk must have a nonzero in-plane 
winding number and hence contains vortices. We denote 
the position of these vortices by k^, and the gaps at the 
vortex cores which determine the Chern number are given 

by 



f* n (k i ) = (g- 1 ) z (e sn (k i )) + h z (k i ). 
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To have a nonzero winding number, the gaps need to have 
different signs at different k^'s. Naively, since (g~ 1 ) z is of 
order U a p, one might expect it to be the dominant term 
and set a large uniform gap. However, in our perturba- 
tion theory, the quasiparticle energies are always close 
to the onsite particle/hole spectrum which are the zeros 
of g _1 (co>), so that as long as the the onsite spectrum 
is nearly spin-symmetric, 25 (g~ 1 ) z (e sn (ki)) ~ 0(h±). 
Therefore, the winding number is determined by the com- 
petition between the terms in Eq. ( [l3| ), and may differ 
than the non-interacting one determined by h z (ki). 

A simple case-To illustrate the relevant physics in 
the simplest setting, we turn off the inter-spin interac- 
tion, setting U-\ = and = U-- = U. We use 

U as the unit of energy and define the dimensionless 
quantities k' = k/y/2mU. The excitation energies are 
£ a (N a ) = UN a — ji a and the zeroth order ground state 
has filling factors N a when N a — 1 < ji a /U < N a . Set- 
ting fi a = U(N a — 1/2)™ the inverse total particle and 
spin propagator is 
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In the cases of equal filling fractions (7V + ,7V_) = 
(TV, TV), (g~ 1 ) z = so that there are no corrections to the 
spin texture. The simplest filling fraction which gives an 
out-of-plane spin texture is (1,0). This filling represents 
a ferromagnetic ground state, which is generally present 
in a bosonic Mott insulator. The quasiparticle disper- 
sions can be expressed implicitly as 
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where e^ sn = (h^ + s\f\^ 3n \)/U, s,n = ±1. Note that 
there is no spin down quasihole band (s,n) = (— , — ) 
because there are no spin down atoms in the unperturbed 
(1,0) ground state. Thus, the ground state momentum 
distribution of spin density is given by Sk = Zk+-fk+-- 

To give an example of hole bands with nonzero Chern 
number, we consider one component of the BHZ Hamil- 
tonian which describes HgTe quantum wells and we 
take parameters such that h = 2t(cosk x + cos/c^), h = 
2t so (sinfc £C , — sin/fy, 0). In the absence of interactions 
these parameters give zero Chern number. However, by 
computing the gaps in Eq. ( 13 ), we find for the filling fac- 
tor (1,2) with t/U ~ .01, and an additional onsite spin 
splitting of S/U ~ .01 that the Chern numbers are ±1, 
so the sum of the Chern numbers is zero for this model. 

An abelian SO coupling-We now consider the spin- 
orbit coupling of Ref. [3 by taking a Hamiltonian that 
at low energy is specified by = k 2 /2m and hk = 
(A x ,r]k x ,0), where m is the effective mass™ f] is the 
spin-orbit coupling and A x is an applied detuning. We 
plot the quasiparticle dispersions in Fig. [lja,b) and Sk 
in Fig. [ljc). The in-plane components of Sk are qualita- 
tively the same as hk, but interactions generate an out- 
of-plane component. Although the non- interacting Berry 
curvature vanishes, h • d x h x d y h = 0, the interacting 
Berry electric field and curvature an d S£ 



„ycx 



respectively, are nonzero, which we plot inFig.[jJd). The 
Chern number is still zero because h has zero in-plane 
winding, and this is verified in our computation because 
is antisymmetric under inversion and hence vanishes 
upon integration. There are still nontrivial physical con- 
sequences due to the quasiparticle anomalous velocity 
= #£ sn (-V0 x z) which gives an edge current in 
the presence of a trapping potential <j). Furthermore, it 
was shown in Ref. [15] for the Fermi liquid that the quasi- 
particle weight is renormalized by a factor Z' ~ — S • V(j). 
Since the only key assumption was a hermitian self en- 
ergy matrix, which is true for our case, we expect this to 
be valid. 

Experimental signatures-There are many experimen- 
tal signatures of the interaction induced Berry curva- 
ture and topological transport that we studied in this 
paper. The ground state particle and spin density can 
be measured using RF spectroscopy.^ and time of flight 
measurements.^ The quasiparticle dispersions and spec- 
tral weights enter into the density-density response func- 
tion which can be measured with Bragg spectroscopy™ 
A nonzero anomalous Hall conductivity causes transverse 
oscillations in collective modes which can be seen by ab- 
sorption imaging of the density. 21 Furthermore, bulk and 
edge quasiparticle transport can be probed directly by 
exciting particle-hole with lasers. 

Conclusion and outlook-We have shown that in a 



spinful bosonic Mott insulator, interactions significantly 
modify the quasiparticles spin-orbit momentum-space 
texture. In our strong-coupling perturbative approxima- 
tion, we have computed the interacting Green function 
for a general spin-orbit coupled hopping Hamiltonian, 
and show that the ground state can in principle support 
an integer Hall conductivity given by the sum of the hole 
Chern numbers. For the optical lattice version of the 
experiment in Ref. [3], we computed quasiparticle dis- 
persions, Berry curvature, and momentum distributions 
of spin density. We hope this will motivate experiments 
to measure these quantities, and to create the spin-orbit 
couplings with in-plane vortex textures which we have 
predicted will exhibit topological transport. 

It seems to be a generic property of two band models, 
which have opposite signs of the Berry charge, that their 
Chern numbers tend to cancel. For future work, we will 
investigate whether it is possible to find a model with 
a nonzero sum of hole Chern numbers, for example, in a 
three-band model. According to the arguments presented 
in Ref. [24], the Hall conductivity for bosons should be 
quantized in even integers. Further research is needed 
to investigate the mutual consistency of these arguments 
with our approximation and the zero temperature Hall 
conductivity that follows from it. 

By choosing the (1,0) filling fraction, we have in a sense 
broken time reversal symmetry "by hand." In the pres- 
ence of inter-spin interactions J7_| , time reversal sym- 
metry will be broken spontaneously because the ground 
state will have magnetic order. 22 This can readily be in- 
cluded in our theory, where Eq. ([9]) still holds but g~ x 
acquires off-diagonal components ™A uniform ferromag- 
netic state gives qualitatively similar results to the case 
studied in this paper. Furthermore, even in the presence 
of magnetization textures, our momentum-dependent re- 
sults will be approximately valid for momenta larger than 
the inverse of the typical length scale of the texture. In 
general, textures will result in interplay between momen- 
tum space and real space berry curvature, which moti- 
vates further study. 

Finally, in this work we have considered only the Mott- 
insulating phase. We leave the analysis of the Mott 
insulator-superfluid transition to future work. We also 
intend to investigate the bulk and edge state transport 
properties of these systems in more detail. 
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